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Abstract 



We investigate a topological inflation model in supergravity. By means 
of numerical simulations, it is confirmed that topological inflation can take 
place in supergravity. We also show that the condition for successful inflation 
depends not only on the vacuum-expectation value (VEV) of inflaton field but 
also on the form of its Kahler potential. In fact, it is found that the required 
VEV of the inflaton tp can be as small as (<p) ~ 1 x Mq, where Mq is the 
gravitational scale. 

PACS numbers: 98.80.Cq,04.65.+e,11.27.+d 
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I. INTRODUCTION 



Superstring theories compactified on (3 + l)-dimensional space-time have many discrete 
symmetries in the low-energy effective Lagrangian [0. A spontaneous breakdown of such 
discrete symmetries creates topological defects, i.e. domain walls, in the early universe [E2|. 
If the vacuum-expectation value (VEV) of a scalar field (p is larger than the gravitational 
scale Mq ~ 2 x 10 18 GeV, the region inside the wall undergoes inflationary expansion and 
eventually becomes the present whole universe ||||] . If the universe is open at the beginning, 
it expands and the spontaneous breakdown of the symmetries always takes place at some 
epoch in the early universe. It has been recently argued that the quantum creation of the 
open universe may take place with appropriate continuation from the Euclidean instanton 
0. Thus, topological inflation is a natural consequence of the dynamics of the system, and it 
does not require any fine-tuning of initial conditions for the beginning universe. Furthermore, 
it does not cause the "graceful exit" problem and the universe becomes homogeneously 
radiation dominated after reheating. 

A simple and interesting model for topological inflation was proposed in the framework 
of supergravity || .[] However, it was not explicitly shown whether topological inflation really 
takes place. In this paper, we perform a numerical analysis on the above model and show 
that topological inflation indeed occurs in a wide range of parameter space. We also show 
that the condition for successful inflation depends not only on the superpotential, which 
determines the vacuum expectation value (VEV) of inflaton ip, but also on the form of its 
Kahler potential. We in fact find that the required VEV can be as small as (ip) ~ 1 x M G , 
which is far below the lower bound of (<p) = r\ cr ~ 1.7M G derived in Ref. ||. 

II. TOPOLOGICAL INFLATION MODEL 

We begin with the topological inflation model proposed in Ref. B, which is based on 
i?-invariant supergravity. The gravitational scale Mq is set to be unity below. In this model 
the superpotential for the inflaton superfield 4>(x, 9) is given by 

W = v 2 X{l-g<p 2 ). (1) 

Here, we have imposed U(1)r x Z 2 symmetry and omitted higher-order terms for simplicity. 
Under the U{1)r we assume 

X{B) -> e- 2ia X(6e ia ), 0(0) -> 0(0e ia ). (2) 

We also assume that the superfield X is even and <fi is odd under the Z 2 . This discrete 
Z 2 symmetry is an essential ingredient for the topological inflation fUJiJ]. In the above 
superpotential (|1]), we always take v 2 and g to be real constants without loss of generality. 
The R- and Z 2 -invariant Kahler potential is given by 

KfaX) = \X\ 2 + \<p\ 2 + hixm 2 + ||V| 4 + ||0| 4 + ■ ■ • , (3) 



1 Other topological inflation models were studied in the superstring inspired models 
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where k±, k 2 and k% are constants of order unity. 

The potential of a scalar component of the superfields X(x, 6) and <p(x, 6) in supergravity 
is given by 

V = eK i(^-^) D Z WD Z *W*-3\W\ 2 } (:, = ,,..Y). (1) 



dzjdz* 



with 



This potential yields an /^-invariant vacuum 

at which the potential energy vanishes. Here, the scalar components of the superfields are 
denoted by the same symbols as the corresponding superfields. If 77 is larger than the critical 
value 7] CT which will be discussed in the next section, the topological inflation occurs. 
For \X\ and \<f>\ <C 1, we approximately rewrite the potential (|J) as 

V ~ v 4 \l - g<p 2 \ 2 + v\l - k^ 2 - k 2 v A \X\ 2 . (7) 

If A)2 ^ —1, X field quickly settles down to the origin and we set X = in our analysis 
taking hi % —1. For g > 0, we can identify the inflaton field tp(x)/y2 with the real part 
of the field 4>(x) since the imaginary part of <p(x) has a positive mass and the real part has 
a negative mass. Because the positive mass of the imaginary part is larger than the size of 
the negative mass, the imaginary part is irrelevant for the inflation dynamics and hence we 
neglect it. Then, we obtain a potential for the inflaton for ip <C 1: 

V(<p) ^v±- ^V, (8) 

where 

K = 2g + k 1 -l. (9) 

The slow- roll condition for the inflaton (p is satisfied for < k < 1 and ^ (p ^ 1.0 The 
Hubble parameter during the inflation is given by H ~ v 2 / v3. The scale factor of the 
universe increases by a factor of e N when the inflaton ip rolls slowly down the potential from 
Pn to 1. The e-fold number iV is given by 

iV~--ln<^v. (10) 

K 

The amplitude of primordial density fluctuations 5p/p due to this inflation is written as 



2 We can always take tp positive since we have the Z 2 symmetry 
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5 -l * 1 ^ „ 2 .o x 10- 5 . (11) 

P 5V37T KlfN 

The normalization is given by the data of anisotropies of the cosmic microwave background 
radiation (CMB) by the COBE satellite [0. Since the e-fold number N corresponding to 



the COBE scale is about 60, which leads to 

w~2.3x 10" 2 v^e-^|jv =60 ~ 1.8 x 1(T 3 - 3.6 x 1(T 4 , (12) 



for 0.02 < k < 0.1. 

The interesting point on the above density fluctuations is that it results in the tilted 
spectrum whose spectrum index n s is given by 

n s ~ 1 - 2k. (13) 

We may expect a possible deviation from the Harrison-Zeldvich scale-invariant spectrum 
n s = 1. Observational constraint on n s is \n s — 1| < 0.2 [|K]], which implies < k < 0.1. 

After inflation ends, the inflaton ip may decay into ordinary particles as discussed in 
Ref. || and the reheating temperature is low enough to avoid overproduction of graviti- 
nos Jl l|] which are thermally produced at the reheating epoch. Recently, nonthermal pro- 



duction at the preheating stage was found to be important in some inflation models |12[. For 
the present model, non-thermal production of gravitinos at the preheating phase is roughly 
estimated as 

(m) „ < 10-14 ( T * ) (14) 

where m v (~ v 2 ), n.3/2, s, and Tr are the mass of the inflaton, the number density of grav- 
itinos, entropy density and reheating temperature, respectively. This is much less than the 
thermal production given by (n 3 / 2 /s)TH ~ 10 _11 (T R /10 10 GeV) and hence we can neglect the 
nonthermal production of gravitinos. 



III. NUMERICAL SIMULATION 

We perform numerical simulations to decide whether topological inflation takes place 
in supergravity and determine the condition for successful topological inflation. For the 
purpose, we follow time evolution of the domain wall and investigate whether it inflates 
or not. Since we consider a planar domain wall whose width is order of the horizon scale, 
we cannot adopt the Friedmann Robert son- Walker metric. Instead, we assume that the 
spacetime has a reflection symmetry of the coordinate x perpendicular to the wall and 
adopt the metric given by 

ds 2 = - dt 2 + A 2 (t, \x\) dx 2 + B 2 (t, \x\) (dy 2 + dz 2 ), (15) 

where A and B correspond to the scale factors in the direction of x and y-z, respectively. If 
the inflation occurs and the proper width of the wall becomes much larger than the horizon 
scale, A and B expand as A ~ B oc e Ht (as shown later) and the universe approaches the 
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de Sitter spacetime. Thus, we examine whether the proper width of the wall becomes much 
larger than the horizon scale. Once it is realized, the universe expands exponentially HQ. 

We adopt the numerical technique developed in Ref. ||. The Einstein-Hilbert action is 
given by 



S 



n 



1 



2 . 2 {d^f-v{ v ) 



(16) 



Variating the above action with respect to the metric g^, we obtain the Einstein equations, 

1 



where Tn U is the energy-momentum tensor, 



-7? =T 

2 /v ± p.vi 



-(p&Y + v{<p) 



(17) 



(18) 



Variation with respect to the scalar field (p gives the equation of motion for the scalar field 



□ 



<P 



dV(<p) 
dip 



(19) 



In order to make it easier to follow time evolution of the system, we choose certain combi- 
nations of Einstein equations, which read 



KL\{2fC — 3/C2) 



2B" 



B 



+ 



2MB' 



A 2 B A 2 B 2 A 3 B 



j2 



01 



■(G\ + G\ + G\ — G 



1 B' 

K>2 + ^"(3^2 ~ 



2 W ' 



B 



,2 



2A 2 B 2 



2\2 



4 AA 2 



V{<p) 



(20) 



(21) 



(22) 



(23) 



where an overdot denotes the time derivative and a dash the spatial derivative. fC{j are the 
extrinsic curvature tensors of constant time hypersurface, given by 



A 



ir 2 — r 3 — __ 
^2 — ^3 — ^1 



(24) 



and K denotes its trace K = K\. The equation of motion for the scalar field <p becomes 



5 



.. „ . if ( A' 2B'\ if dV(ip) „ 

v _^_^_(__ + _j|_ + _ip = 0. (25) 

We set the initial condition for numerical simulations. First, we consider an initial 
configuration of a domain wall. For the convenience of numerical calculations, we take only 
the region \x\/S < 2 [] where S(— is the width of the domain wall and 77 is the VEV 

of if. We impose the free boundary condition, that is, ip' = A' = B' = at the boundaries 
x/5 = —2 and x/5 = 2. Then, we adopt the following initial configuration for the domain 
wall so that the gradient of the field disappears at the boundaries, 



ip(t = 0, x) 



7/ 



x 5/8 x\ 3 3 / 8 x N v 
6 ~ 4 KTEl) + 8 V15 5 



x 15\ 
15 x 



with ip(t = 0, — x) = —ip(t = 0, x) for —2 < x/5 < 0. This is a deformed version of the static 
domain wall solution in a flat spacetime, </9fl at = i]t&nh(x/5). The function = 0,x) is 
decided so as to satisfy the following three conditions: (1) it is a fifth-order odd polynomial 
function of x, (2) the first term coincides with that of the expansion of (fast, (3) it is smooth 
at x/S — 15/8, that is, ip' = ip" = at x/5 = 15/8. Also, <p is set to be 0. 

Next, on the initial hypersurface, we determine A, B, /C|, K, so as to satisfy the Hamil- 
tonian constraint (|20|) and the momentum constraint (^T]) . We have freedom for the initial 
hypersurface to have homogeneous and isotropic curvature, which automatically satisfies the 
momentum constraint fl2~lD . This choice leads to 

— = K,\ = K\ = "negative" const, (27) 

where "negative" implies that the universe is in an expanding phase. Furthermore, we 
can take the conformally flat spatial gauge, A = B, on the initial hypersurface and set 
A = B = 1 at 1 = 0. Finally, we determine the negative value of K,. Since we adopt 
the reflection symmetry of the coordinate x and the free boundary condition, the condition 
A' = B' = at x = and x = ±25 must be satisfied. K, is determined so that the 
Hamiltonian constraint ( p0[) satisfies the above conditions. 

Now the initial settings are completed and hence we have only to follow the time evolution 
of five variables, A, B, JC, K%, and <p. Note that we have introduced five variables, A, B, JC, K\, 
and ip though only three variables are independent. This is partly because the second order 
differential equations have been reduced to the first order differential equations. Moreover, 
as for time evolution of K\, we use Eq . fl23|) only at x = and acquire the value at x 7^ by 
integrating Eq.Q2"T|) in the direction of x in order to avoid numerical instability. 

When inflation takes place, A and B grow exponentially so that the proper distance 
from the domain wall core (x — 0) also increases exponentially. In order to see whether this 
happens or not, we follow time evolution of the width of the wall for a given potential V(ip). 



3 We have confirmed that the results do not change even if we take wider ranges of the direction 
x (e.g. \x\/5 = 3). 
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To fix the potential V(p), we first consider the Kahler potential with only terms up to 
the fourth order, 

K(<j>, X) = |X| 2 + |0| 2 + hlxm* + ||X| 4 + (28) 
Then, the Lagrangian density is given by 
£(0) = A,n(0) - V{cj>) 

exp(|0| 2 + ^|0 4 f) 

= -(i + m 2 )d,W -v*\i- # 2 | 2 (29) 

Here we have set X = 0. Identifying the inflaton field (p(x)/\/2 with the real part of the 
field <p(x), the Lagrangian density becomes 

C(<p) = Ckin(p) ~ V(tp) 



--A^- 2 ^)^y-^-^) V + v • (30) 



with the VEV r\ — J 2/g. In the present model we have four free parameters k\, k 3 and g. 
However, ^(^ —1) only works as a stabilizer of the X field as explained before and it is not 
important for the dynamics of topological inflation itself. Once the X field is stabilized at 
X = 0, the potential V(ip) does not depend on fc 2 . k 3 is almost irrelevant for the dynamics 
of (f) field and only changes its VEV due to the redefinition of p with a canonical kinetic 
term. Then, we set k$ = first and later consider the case of nonzero k^. Thus, we have 



only two relevant parameters, fci and g. The potential V((p) has a pole at <p = y2/\ki\ for 
k\ < 0. But, we are only interested in the dynamics of ip up to the VEV 77 so that there is 
no problem if \k\\ < g for k\ < 0. 

We introduce dimensionless quantities, (p = ip/M G , x = xH(x — 0), t — tH(x = 0), K.^ = 
K,ij/H{x = 0), and 5 = 5H(x = 0) = l/(y/5g) where H(x = 0) = v 2 /V3. As the first step, 
we consider the simplest case of k\ = and g = 0.5 (i.e., 77 = 2.0), which leads to the spectral 
index n — 1. Time evolution of the domain wall is depicted in Fig. [1|. The vertical axis 
represents the value of the scalar field. The horizontal axis represents the proper distance 
from the domain wall core. As time elapses, the domain wall (the region for ip < 0.8) 
expands and topological inflation really takes place. As shown in Figs. || and |3], once the 
domain wall expands enough, the scale factors A and B increase at the same rate, that is, 
— K\ ~ —K>2 ~ H(x = 0) inside the wall. Thus, the universe expands exponentially and 
approaches the de Sitter spacetime. 

We consider the dependence on k for the VEV fixed, rj = 2.0 (g = 0.5). The result is 
depicted in Fig. |j. For large k, the domain wall cannot inflate enough, k controls the mass 
scale of the scalar field p near the origin. As k becomes larger, the scalar field tp rolls down 
faster so that only the small region of the original domain wall earns the vacuum energy and 
cannot overcome the gradient energy. 

The dependence on v is studied, which determines the energy scale of the domain wall. 
The result with the same parameters in Fig. [| except for v = 3.6 x 10~ 4 is shown in Fig. 
|5|. The results are quite the same and have no dependence on v. This can be interpreted as 
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follows: First, the dependence on v only appears through the Hubble parameter during the 
inflation given by H ~ v 2 /y/3. But, since the width of the domain wall 5 ~ l/(y/3g~H), the 
rough criterion 5 > H~ l becomes independent of H, that is, v. Also, rewriting the potential 
as V((p) ~ \((p 2 — rj 2 ) with A = (v/rj) 4 , the independence of v is equivalent to that of A. The 
independence of v is also found in Ref. 0. 

We now consider the cases of the nonzero k 3 . Since the kinetic term Chin is not canonical 
in these cases, we define the scalar field $ with the canonical kinetic term as 



1 / 1 arcsinh ( J^(p ) 



MM, 



(31) 



1 / i arcsm l y H^tp 
-<p\l - -\k 3 \tp 2 + ^ (h<0). 



The results for the same parameters in Fig. [l| except for k 3 = ±0.3 are depicted in Fig. 
^|. The positive k 3 encourages the occurrence of topological inflation while the negative k 3 
discourages. This is because the VEV of $ is larger than rj for the positive k 3 while smaller 
for the negative k 3 . 

Finally, we discuss the criterion for the VEV of (p for successful topological inflation 
within < k ^ 0.1. In previous analyses, we search for only the parameter region satisfying 
| fci| < g for negative k\ because of the appearance of a pole of the potential. The condition 
\k\\ < g leads to g < 1 + k so that rj > 1.41(1.35) for k = 0.0(0.1), where we have confirmed 
the occurrence of topological inflation. In order to obtain the lower limit of rj, we add to 
the Kahler potential the sixth order terms, 

AK = hlXl 2 ^ + l 2 \Xm 2 + ||V| 6 + ^|0| 6 , (32) 

where there are four parameters but only l\ is relevant. If we take l\ satisfying l\ > g(g — 1), 
the pole smaller than the VEV does not appear in the whole k\ — g parameter space. The 
result for smaller rj with k = 0.0 is depicted in Fig. [7|. We find the critical value of the 
breaking scale, i] cr ~ 0.95 (1.00) for k = 0.0 (0.1), which corresponds to (4>) ~ l/\/2. 



IV. CONCLUSIONS AND DISCUSSIONS 

We have studied a topological inflation in supergravity. First, we have shown that 
topological inflation really takes place in supergravity. Also, the criterion of successful 
topological inflation depends not only on the breaking scale of the discrete symmetry but 
also on the mass of the inflaton near the origin. This is because the inflaton rolls down 
rapidly from the origin if its mass is large. For a very flat case favored by the observation of 
the spectral index, n s ~ 1 — 0.8 (i.e., < k < 0.1), we have found that the critical breaking 
scale f] cr becomes as small as Mq, which is smaller than the critical value, rj cr = 1.7 Mq 
observed in Ref. ||. Finally we have discussed the primordial spectrum produced by the 
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topological inflation. In general, the topological inflation predicts the tilted spectrum n s < 1 
depending on 

The present topological inflation model is free from the thermal and nonthermal over- 
production of gravitinos since the reheating temperature can be as low as 10 8 GeV. Further- 



more, as pointed out in Ref. ||13|] , this model is consistent with a leptogenesis scenario in 
which heavy Majorana neutrinos are produced in the inflaton decay and successive decays of 
the Majorana neutrinos result in lepton asymmetry enough to explain the observed baryon 
asymmetry in the present universe. 
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4 It is possible to produce more exotic spectrum including blue one. This is due to the exponential 
blow of the potential, which is significant for the region ip > Mq and makes the potential more 
complex than the simple double-well potential. 
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FIG. 1. Time evolution of the domain wall in the case of k\ = and g = 0.5 (rj = 2.0). The 
vertical axis represents the value of the scalar field. The horizontal axis represents the proper 
distance from the domain wall core. Note that the proper width of the wall becomes much larger 
than the Horizon scale. As time elapses, the domain wall expands and topological inflation takes 
place. Here we set v = 1.8 x 10 -3 . But the result does not depend on the energy scale v as shown 
later. 
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FIG. 2. The expansion rates of the scale factors, A and B, that is, — /C{ = A/ A and — = B/B 
are shown for the case of Fig. [l]. As the universe expands enough, — K,\ and — /C| approach the 
same value, the Hubble parameter = 0), inside the wall. 
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FIG. 3. Time evolution of — /C} and — /C| at the origin x = is shown for the case of Fig. 
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FIG. 4. Time evolution of the domain wall in the cases of different k for the fixed VEV, 
77 = 2.0 (g = 0.5), with v = 1.8 X 10~ 3 . 
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FIG. 5. The results with v = 3.6 x 10" 4 . The energy scale v is different from that in Fig 
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FIG. 6. The results for the same parameters in Fig. [T] except k% = ±0.3. 
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FIG. 7. The results with rj ~ r\ cr for k = 0.0. They show r\ cr ~ 0.95 for k = 0.0 
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